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ABSTRACT  Background:  Repetitive  hypobaric  exposure  in  humans  induces  subcortical  white  matter  change, 
observable  on  magnetic  resonance  imaging  (MRI)  and  associated  with  cognitive  impairment.  Similar  findings  occur  in 
traumatic  brain  injury  (TBI).  We  are  developing  a  swine  MRI-driven  model  to  understand  the  pathophysiology  and  to 
develop  treatment  interventions.  Methods:  Five  miniature  pigs  ( Sus  scrofa  domestica)  were  repetitively  exposed  to 
nonhypoxic  hypobaria  (30,000  feet/FI02  100%/transcutaneous  P02  >90%)  while  under  general  anesthesia.  Three  pigs 
served  as  controls.  Pre-exposure  and  postexposure  MRIs  were  obtained  that  included  structural  sequences,  dynamic 
contrast  perfusion,  and  diffusion  tensor  quantification.  Statistical  comparison  of  individual  subject  and  group  change 
was  performed  utilizing  a  two-tailed  t  test.  Findings:  No  structural  imaging  change  was  noted  on  T2-weighted  or 
three-dimensional  fluid-attenuated  inversion  recovery  imaging  between  MRI  1  and  MRI  2.  No  absolute  difference  in 
dynamic  contrast  perfusion  was  observed.  A  trend  ( p  =  0.084)  toward  increase  in  interstitial  extra- axonal  fluid  was 
noted.  When  individual  subjects  were  examined,  this  trend  toward  increased  extra-axonal  fluid  paralleled  a  decrease  in 
contrast  perfusion  rate.  Discussion/Impact/Recommendations:  This  study  demonstrates  high  reproducibility  of  quantita¬ 
tive  noninvasive  MRI,  suggesting  MRI  is  an  appropriate  assessment  tool  for  TBI  and  hypobaric-induced  injury 
research  in  swine.  The  lack  of  fluid-attenuated  inversion  recovery  change  may  be  multifactorial  and  requires  further 
investigation.  A  trend  toward  increased  extra-axonal  water  content  that  negatively  correlates  with  dynamic  contrast  per¬ 
fusion  implies  generalized  axonal  injury  was  induced.  This  study  suggests  this  is  a  potential  model  for  hypobaric- 
induced  injury  as  well  as  potentially  other  axonal  injuries  such  as  TBI  in  which  similar  subcortical  white  matter  change 
occurs.  Further  development  of  this  model  is  necessary. 


INTRODUCTION 

Increased  subcortical  white  matter  hyperintensity  (WMH) 
burden  on  magnetic  resonance  imaging  (MRI)  is  associated 
with  nonhypoxic  hypobaric  occupational  exposure  in  U-2 
pilots  and  in  altitude  chamber  inside  safety  monitors.1-3 
Additionally,  divers  have  increased  WMH  burden.4  Extreme 
mountain  climbers  have  white  matter  (WM)  volume  change5 
and  diffusion  tensor  imaging  (DTI)  changes6,7  and  U-2  pilots 
have  DTI  change  after  excluding  contribution  from  WMH.8 
Recently,  we  have  demonstrated  that  a  single  hypobaric 
exposure  in  humans  does  not  induce  acute  WM  hyperinten¬ 
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sities  on  fluid-attenuated  inversion  recovery  (FLAIR)  but  is 
associated  with  other  transient  MRI  changes  (P.  Sherman, 
S.  McGuire,  personal  communication;  August  2,  2016).  All 
of  these  populations  experience  decompressive  stress  and 
are  at  risk  of  decompression  sickness.  Reduced  executive 
cognitive  performance  has  been  linked  to  WMH  burden9,10 
and  disruption  of  axonal  integrity  as  measured  by  DTI.11,12 
The  traditional  neuropathological  model  posits  gaseous 
embolic  damage  to  the  brain  as  the  main  mechanism  respon¬ 
sible  for  hypobaric-induced  brain  injury,  although  the 
possibilities  of  microemboli  and/or  inflammatory  and/or 
immunological  etiologies  inciting  a  more  diffuse  process 
have  been  suggested.13-16  Given  the  presumed  permanence 
of  the  subcortical  WM  brain  injury  and  reduced  neuro- 
cognitive  performance  associated  with  hypobaric  exposure, 
an  animal  model  is  necessary  to  understand  the  pathophysi¬ 
ology  and  to  develop  prevention  and  treatment  interventions. 

We  are  developing  a  translational  neuroimaging  animal 
model  to  understand  the  physiology  of  diffuse  brain  injury 
following  nonhypoxic  hypobaric  exposure.  We  chose  the 
miniature  pig  ( Sus  scrofa  domestica)  as  the  animal  model 
because  of  the  close  resemblance  of  the  pig  gyrencephalic 
brain  to  the  human  brain  in  both  white-to-grey  matter  ratio 
and  in  cerebral  circulation.  Likewise,  the  same  MRI  scanner 
and  sequences  can  be  used  in  both  humans  and  miniature 
pigs,  providing  a  translational  component  to  this  research. 
Although  large  animals  with  a  greater  degree  of  WM  would 
have  been  ideal,  the  size  of  the  hypobaric  chamber  and  the 
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desire  to  employ  the  same  MRI  sequences  utilized  in  our 
human  studies3  restricted  the  animal  model  to  a  weight  of 
20  to  25  kg. 

We  hypothesized  that  the  MRI  changes  in  this  model 
would  parallel  the  changes  reported  in  humans.  We  studied 
diffusion  changes  in  cerebral  WM  using  the  T2-weighted 
three-dimensional  (3-D)  FLAIR,  dynamic  susceptibility  con¬ 
trast  (DSC),  and  advanced  diffusion-weighted  signal  with  a 
longitudinal  imaging  approach.  DSC  utilizes  a  contrast  agent 
to  assess  vascular  integrity.  Advanced  diffusion- weighted  sig¬ 
nal  (multi-Z?-value  and  permeability-diffusivity  [PD])  assesses 
a  two-compartment  restricted  and  unrestricted  water  model  to 
assess  changes  in  interstitial  water  content.  As  the  hypobaric 
injuries  in  humans  are  mainly  observed  in  the  cerebral  WM, 
the  amount  of  WM  in  the  brain  of  the  miniature  pig  compared 
to  rodent  models  would  increase  the  likelihood  of  finding 
WMH  on  MRI  and  neuropathological  examination. 

METHODS 

Animal  Subjects 

All  animal  protocols  were  reviewed  and  approved  by  the 
Institutional  Animal  Care  and  Use  Committee,  59th  Medical 
Wing,  Joint  Base  San  Antonio-Lackland  AFB,  Texas.  Sev¬ 
enteen  miniature  female  pigs  (14  study  subjects;  3  control 
subjects)  weighing  19  to  24  kg  were  used.  Six  subject  ani¬ 
mals  were  lost  to  severe  pulmonary  injury  from  ventilator 
barotrauma  and/or  severe  pulmonary  decompression  sickness 
with  subsequent  overwhelming  sepsis  before  MRI  2,  and 
two  sustained  cardiac  arrest  during  altitude  exposure  and 
one  experienced  a  fatal  anaphylactic  allergic  reaction  to 
propofol,  leaving  eight  animals  available  for  baseline  and 
postexposure  MRI  assessment  (five  subjects  and  three  con¬ 
trols;  Table  I).  Only  female  pigs  were  used,  as  bladder  cath¬ 
eterization  was  necessary  for  fluid  monitoring  during  the 
prolonged  periods  of  anesthesia  for  the  chamber  exposures 
and  for  the  MRIs.  All  animals  underwent  the  same  proce¬ 
dures,  including  placement  in  the  barometric  chamber  for 
8  hours  while  under  general  anesthesia,  the  only  difference 
being  hypobaric  exposure.  The  three  control  animals  were 
maintained  at  sea  level  on  room  air  in  the  chamber,  whereas 
the  subject  animals  were  taken  to  9,144  m  (30,000  feet)  on 


100%  oxygen.  All  imaging  was  performed  at  the  Wilford 
Hall  Ambulatory  Surgical  Center,  59th  Medical  Wing.  Ani¬ 
mals  were  housed  in  the  center’s  animal  facilities. 

Study  Protocol 

Anesthesia  during  imaging  was  maintained  using  isoflurane 
via  an  MRI-compatible  anesthesia  machine.  Since  gaseous 
anesthesia  was  unreliable  at  altitude,  anesthesia  during 
active/sham  hypobaric  exposure  was  maintained  using  peri¬ 
odic  boluses  of  propofol/ketamine  adjusted  to  maintain 
stable  physiological  parameters  and  anesthesia.  On  study 
day  1,  baseline  imaging  was  performed.  Exposure  episodes 
began  on  study  day  3  with  three  subject  animals  exposed  six 
times  to  9,144  meters  (ascent/descent  time  15  minutes)  over 
12  days,  one  exposed  five  times  and  one  exposed  four  times. 
Three  control  animals  underwent  six  sham  exposures  at  sea 
level.  Following  intubation,  all  animals  received  1  hour  of 
100%  oxygen  prebreathe  before  active  or  sham  exposure. 
Animals  taken  to  altitude  remained  on  100%  oxygen  and  the 
sham  animals  were  placed  on  room  air  at  sea  level.  As  venti¬ 
lator  adjustments  were  not  possible  once  a  subject  animal 
was  in  the  chamber,  tidal  volumes  were  preset  on  the  basis 
of  weight  and  volume  calculations  derived  from  Boyle’s 
law.  Throughout  exposure,  transcutaneous  pulse  oxygen 
>90  mmHg  and  endotracheal  tube  carbon  dioxide  35  to 
45  mmHg  were  maintained  by  drug  infusion  adjustment. 
Chamber  temperature  was  maintained  between  68  and  70°F. 
Fluid  input  was  adjusted  to  match  urine  output.  On  study 
day  15  (2  days  after  final  exposure),  a  second  MRI  was 
performed.  Subsequently,  all  animals  were  euthanized  with 
sodium  pentobarbital.  Necropsy  including  neuropathological 
examination  of  neural  tissue  was  performed. 

Magnetic  Resonance  Imaging 

All  imaging  data  were  collected  at  the  Wilford  Hall  Ambulatory 
Surgical  Center,  59th  Medical  Wing,  Joint  Base  San  Antonio- 
Lackland,  Texas,  using  a  Siemens  3T  Verio  scanner  (Siemens 
Medical  Solutions  USA,  Inc.,  Malvern,  Pennsylvania)  equipped 
with  a  15 -channel  phased  array  knee  coil  applied  as  a  head 
coil  for  pigs  in  the  prone  position.  The  protocol  consisted  of 
high-resolution  T2-weighted  FLAIR  sequence,  multi-Z?-value 
diffusion- weighted  imaging  protocol,  and  DSC. 


TABLE  I.  Subject  Specifics 


Subject 

Weight/Sex  (kg/F) 

Flights 

Cause  of  Death 

Control  1 

19 

6 

Euthanized  after  MRI  2 

Control  2 

19 

6 

Euthanized  after  MRI  2 

Control  3 

22.5 

6 

Euthanized  after  MRI  2 

Subject  1 

17 

6 

Barotrauma-induced  death  3  days  after  MRI  2 

Subject  2 

20 

4 

Euthanized  after  MRI  3 

Subject  3 

22 

5 

Euthanized  after  MRI  3 

Subject  4 

24 

6 

Euthanized  after  MRI  3 

Subject  5 

19 

6 

Fatal  sepsis  resulting  from  severe  sinusitis  6  days  after  MRI  2 
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T2-Weighted  FLAIR  Imaging  and  Analysis 

Three-dimensional  FLAIR  provides  structural  imaging  of  the 
brain  with  the  suppression  of  cerebrospinal  fluid  signal.17 
Imaging  was  performed  in  a  similar  manner  to  previously 
reported  human  studies3  except  with  an  isotropic  resolution  of 
0.30  x  0.30  x  0.30  mm.  A  3-D  turbo  spin-echo  sequence  used 
the  following  parameters:  TR  =  4,500  milliseconds,  TE  = 
311  milliseconds,  flip  angles  =  90/180,  echo  train  length  of 
224,  and  nonselective  inversion  pulse  with  TI  =  1,800  milli¬ 
seconds.3  Mango  software  (http://ric.uthscsa.edu/mango)  was 
used  to  display  the  FLAIR  images,  permitting  semiautomatic 
delineation  of  WMH  lesions  by  an  experienced  neuroanatomist 
as  previously  described.2 

Multi-b-Value  Imaging  Protocol  and  PD  Model  Calculations 
This  protocol  was  developed  on  the  basis  of  q-space  in  vivo 
mapping  of  water  diffusion.18-22  This  protocol  consisted  of 
15  shells  of  b  values  (b  =  250,  500,  600,  700,  800,  900, 
1,000,  1,250,  1,500,  1,750,  2,000,  2,500,  3,000,  3,500,  and 
3,800  s/mm2;  diffusion  gradient  duration  =  47  milliseconds, 
separation^  54  milliseconds).  Thirty-five  isotropically  dis¬ 
tributed  diffusion-weighted  directions  were  collected  per 
shell,  including  16  b  =  0  images.  The  highest  b  value  ( b  = 
3,800  s/mm2)  was  chosen  because  the  signal-to-noise  ratio 
(SNR)  for  the  corpus  callosum  (SNR  =  6.1  ±  0.7)  during 
protocol  development  met  the  empirical  requirement  of 
SNR  =  5.0.  The  b  values  and  the  number  of  directions  per 
shell  were  chosen  for  improved  fit  of  the  biexponential 
model.23  The  data  were  collected  using  a  single-shot,  echo- 
planar,  single  refocusing  spin-echo,  T2-weighted  sequence 
(TE/TR  =  120/1,500  millisecond  with  the  FOV  =  100  mm) 
with  a  spatial  resolution  of  1  x  1  x  3  mm  and  12  slices  pre¬ 
scribed  in  axial  orientation  to  sample  the  entire  brain.  The 
scan  time  was  about  10  minutes. 

The  WM  mask  was  derived  per  subject  on  the  basis  of 
the  contrast  in  fractional  anisotropy  (FA)  between  WM, 
GM,  and  CSF.  Voxel-wise  FA,  DA,  and  DR  images  were 
created  using  Camino  software  (http://cmic.cs.ucl.ac.uk/ 
camino).24  Segmentation  was  semiautomatic  for  each  ani¬ 
mal,  using  an  intensity  histogram  approach  and  manual 
editing  in  the  Mango  software. 

The  PD  model  was  presented  elsewhere  25  using  parameters 
derived  from  the  following  biexponential  two-compartment 
diffusion  model  fit: 

Here  S(Z?)  is  the  average  diffusion- weighted  signal  for  a 
given  b  value,  averaged  across  all  directions.  This  model 
is  focused  on  calculation  of  Mu,  which  is  the  compart- 
mental  fraction  of  the  signal  that  comes  from  interstitial 
(unrestricted)  water,  and  (1  -  Mu)  is  the  fraction  of  the 
compartment  with  restricted  diffusion  (Mr).  This  model 
assumes  that  the  diffusion  signal  is  produced  by  two 
quasi-pools  of  anisotropically  diffusing  water.  Du  is  the 
mean  unrestricted  diffusivity  of  the  water  molecules  that 
are  away  from  the  axonal  membranes.  Dr  is  the  mean 


restricted  diffusivity  of  the  water  molecules  that  are  within 
the  restricted  ion  channels  or  close  to  axonal  membranes. 
The  parameters  are  the  compartmental  fraction  Mu  and  the 
PD  index  (PDI),  which  is  the  ratio  of  Dr  and  Du.  An 
increase  in  Mu  would  correspond  to  the  rise  in  the  fraction 
of  interstitial  water  within  WM  voxels  that  can  be  inter¬ 
preted  as  unrestricted  water  molecules  that  are  away  from 
the  axonal  membrane. 


DSC  Perfusion  Imaging  and  Analysis 

DSC  perfusion  imaging  is  a  sensitive  approach  for  the  evalu¬ 
ation  of  embolic  cerebral  damage  and  change  in  permeabil¬ 
ity  of  the  blood  brain  barrier  (BBB). 26,27  DSC  utilizes  an 
exogenous  contrast  imaging  agent  that  historically  has  been 
used  to  study  changes  in  cerebral  perfusion  and  disruptions 
of  the  BBB  caused  by  trauma,  ischemia,  cerebral  inflamma¬ 
tion,  multiple  sclerosis,  neoplasms,  and  other  conditions.28,29 
Under  a  normal  physiological  state,  contrast  agents  used 
with  DSC,  such  as  diethylenetriamine  pentaacetic  acid  bis- 
methylamide  chelate  of  gadolinium  (Gd-DTPA;  Magnevist; 
molecular  weight  573.3  Da)  (Berlex  Laboratories,  Wayne, 
New  Jersey),  cannot  cross  the  BBB.  If  the  BBB  is  disrupted, 
the  contrast  agent  enters  surrounding  tissues  and  becomes 
temporarily  entrapped  in  the  cellular  matrix.  This  entrapment 
is  detected  by  alterations  in  the  washout  rate  and  reduced 
intensity  in  the  affected  tissues;  specifically,  on  a  suscepti¬ 
bility  (T2*)  sensitive  imaging  sequence,  Gd-DTPA  alters 
the  magnetic  relaxation  properties  of  the  blood  and  therefore 
elongates  washout  curves,  preventing  full  recovery  of  imag¬ 
ing  intensity.30  The  signal  intensity  is  reduced  in  linear 
proportion  to  the  concentration  of  the  contrast  agent  and, 
therefore,  produces  regional  reductions  in  signal  proportional 
to  the  perfusion  rates. 

A  10-mL  bolus  (46.901  mg/mL;  mean  23.3  mg  Gd-DTPA/ 
kg  of  body  weight;  range  19.5-27.6  mg)  of  Gd-DTPA  was 
administered  to  the  pig  through  a  port-a-catheter  inserted 
1  week  before  experiment.  The  baseline  and  postexposure 
imaging  protocols  included  a  single  DSC  study  utilizing  a 
gradient  echo-planar-imaging  perfusion  sequence.  Parame¬ 
ters  were  adjusted  to  obtain  high  temporal  (TR  =  2)  and 
spatial  (1.6  x  1.6  x  2.0  mm)  resolution,  acquiring  23  slices 
for  full  brain  coverage  with  no  gaps.  Imaging  slices  were 
prescribed  axially. 

Image  preprocessing  was  performed  using  Functional 
Magnetic  Resonance  Imaging  Brain  Software  Library  (http:// 
fsl.fmrib.ox.ac.uk/fsl/fslwiki/)  and  consisted  of  motion  cor¬ 
rection  and  automated  removal  of  nonbrain  tissue.  Following 
this,  the  average  image  was  loaded  in  Mango  software 
where  the  remaining  nonbrain  tissue  was  removed  via  man¬ 
ual  detailing. 

Mango  was  used  to  calculate  the  average  brain  intensity 
measurements  for  each  of  the  time  points  of  DSC  sequence. 
Value  normalization  was  performed  by  averaging  all  intensi¬ 
ties  before  injection  (mean  data  points  averaged  64,  range 
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51-65)  to  produce  an  average  baseline.  The  measured  inten¬ 
sities  were  then  normalized  using  the  equation: 

A  7  Measured  Intensity  —  Min  Intensity 

Normalized  =  — - - - - - 

Avg Initial  Baseline  —  Min  Intensity 

where  “Normalized  Intensity”  is  the  value  corresponding  to 
the  normalized  0-1  scale,  “Measured  Intensity”  is  the  origi¬ 
nal  intensity  value,  “Avg  Initial  Baseline”  is  the  average  of 
the  first  51  to  65  points  before  injection  of  the  Gd-DTPA, 
and  “Min  Intensity”  is  the  minimum  intensity  following  the 
injection  of  Gd-DTPA.  Normalized  signal  intensity  trends 
were  shifted  to  have  the  minimum  value  at  time  point  72 
(140  seconds)  to  accommodate  for  potential  variability  in 
the  time  of  injection. 

The  recovery  of  signal  intensity  corresponding  to  the 
washout  of  contrast  agent  from  the  brain  was  fit  using  a 
natural  logarithmic  function:  Intensity  =  A  x  ln(/)  +  B, 
where  t  is  the  time  and  A  and  B  are  fit  constants.  For  this 
equation,  t  =  0  corresponds  to  the  minimum  intensity  at 
point  72  (140  seconds).  The  washout  rate  was  derived  from 
the  average  of  the  first  derivative  values  of  the  functional 
model  for  the  initial  washout,  using  the  equation: 

dl  _  A 
dt  t 

where  dlldt  is  the  first  derivative,  A  is  the  fit  constant  from 
the  original  equation,  and  t  is  time.  The  postcontrast  inten¬ 
sity  was  the  average  of  the  final  function  values.  The  wash¬ 
out  rate  and  final  value  were  used  to  evaluate  changes 
between  the  first  and  second  imaging  of  each  subject.  The 
intensities  for  the  sham  and  active  animals  were  averaged 
respectively,  and  the  same  procedure  was  applied. 


Statistical  Analysis 

Individual  subject  baseline  and  postexposure  results  were 
compared  using  2-tailed  paired  t  tests  to  evaluate  changes 
between  the  first  and  second  imaging  of  each  subject.  Addi¬ 
tionally,  a  2-tailed  normal  t  test  was  used  for  group  compari¬ 
son  of  sham  and  active  subject  results.  We  used  chi-squared 
and  Spearman’s  correlation  factor  for  assessment  of  validity 
of  observed  versus  derived  equation  equivalency.  We  con¬ 
sidered  p  <  0.05  to  be  significant. 

RESULTS 

There  was  no  group  difference  in  the  body  weight  between 
sham  and  active  animals  (p  =  0.89).  No  WMH  lesions  or 
vasogenic  edema  changes  were  detected  on  the  baseline  or 
follow-up  T2- weighted  or  3-D  FLAIR  images  (Fig.  1). 

There  were  no  significant  changes  in  the  Mu  in  sham  or 
PDI  coefficients  in  active  and  sham  animals,  although  there 
was  a  trend  toward  an  increase  in  Mu  between  the  baseline 
and  follow-up  in  active  subjects  (Table  II;  p  =  0.084). 

A  tight  fit  of  the  observed  averaged  DSC  data  to  the 
exponential  function  was  present  (Fig.  2).  There  were  no 
significant  differences  in  the  fit  parameters  between  baseline 
and  follow-up  imaging  sessions  in  both  sham  and  active  ani¬ 
mals  (Table  III).  Likewise,  there  were  no  significant  differ¬ 
ences  between  baseline  and  follow-up  measurements  of  the 
washout  rate  and  the  ratio  of  final/baseline  intensity  parame¬ 
ters  for  the  sham  or  active  animals  (Table  III).  Finally,  com¬ 
paring  the  washout  rate  and  the  ratio  final/baseline  intensity 
for  sham  to  active  animals,  no  significant  differences  were 
found  between  the  baseline  and  follow-up  imaging  between 
groups  (all  p  >  0.05). 

Correlation  of  DSC  to  Mu  suggested  an  inverse  relationship 
between  washout  rate  and  Mu,  with  increased  unrestricted  inter¬ 
stitial  water  associated  with  decreased  Gd-DTPA  elimination 


FIGURE  1.  Pre-  and  postexposure  3-D  FLAIR  image.  Coronal,  sagittal,  and  axial  representative  images  demonstrate  clear  separation  of  white  and  gray 
matter.  No  WM  hyperintense  lesions  were  apparent  in  the  subcortical  and  periventricular  regions. 
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TABLE  II.  PD-Model  Analysis  of  Diffusion-Weighted  Data 


Item 

MRI  1  (Pre-exposure) 

MRI  2  (Postexposure) 

p  Value  Pre/Post 

Mu 

PDI 

Mu 

PDI 

Mu 

PDI 

Sham  Average 

0.5874 

0.0898 

0.5937 

0.0935 

0.128 

0.136 

SD 

0.0047 

0.0013 

0.0015 

0.0011 

Cl 

0.0053 

0.0014 

0.0017 

0.0012 

Study  Average 

0.5897 

0.0865 

0.5985 

0.1076 

0.084 

0.378 

SD 

0.0054 

0.0083 

0.0102 

0.0460 

Cl 

0.0048 

0.0073 

0.0090 

0.0403 

Cl,  95%  confidence  interval;  SD,  standard  deviation. 


(sham  animals:  y  =  -51.538x  +  56.699;  R 2  =  0.8979;  study 
subjects  y  =  -2.7973x  +  5.4942;  R 2  =  0.0194). 

Neuropathological  examination  demonstrated  neuropil 
clear  spaces,  gliotic  changes,  subarachnoid  and  WM  micro¬ 
hemorrhages,  thoracic  spinal  cord  perivascular  micro¬ 
hemorrhage,  arachnoidal  inflammation,  focal  red  Purkinje 
cells,  focal  collection  of  oligodendrocytes,  and  focal  bacte¬ 
rial  inflammation. 

DISCUSSION 

Our  study  demonstrated  an  excellent  stability  of  the  quanti¬ 
tative  MRI  measurements  in  both  active  and  sham  animals 
and  the  feasibility  of  this  approach  in  a  miniature  pig  model 
of  hypobaric  brain  injury.  The  hypobaric  exposure  demon¬ 
strated  no  observable  WMH  changes  in  the  active  animals. 
WMH  are  nonspecific31-35  but  have  been  associated  in  two 
human  populations  exposed  to  hypobaria  after  controlling 
for  other  potential  etiologies.3  The  failure  to  demonstrate 
formation  of  WMH  in  this  study  may  be  secondary  to 
insufficient  exposure  stimulus,  to  lack  of  sufficient  follow¬ 


up  for  WMH  development,  or  to  inherent  resistance  of  mini¬ 
ature  swine  to  developing  cerebral  WMH  lesions.  The  small 
size  of  this  study  further  limits  the  demonstration  of  signifi¬ 
cance;  human  studies  on  hypobaric-associated  WMH  exam¬ 
ined  83  to  106  subjects  per  group.3  Postmortem  examination 
in  this  study  noted  a  number  of  nonspecific  changes,  but 
how  these  potentially  relate  to  WMH  or  directly  to  hypobaric 
exposure  is  unclear.  Pathological  correlation  to  WMH  in 
humans  has  occurred36  but  to  our  knowledge  this  has  not 
occurred  in  swine. 

No  statistically  different  changes  were  noted  in  DSC.  A 
trend  was  noted  in  the  fraction  of  unrestricted  interstitial 
(extra-axonal)  water  (Mu)  on  the  analysis  of  the  diffusion- 
weighted  images.  Furthermore,  there  was  an  inverse  corre¬ 
lation  of  Mu  with  Gd-DTPA  washout,  suggesting  injury 
had  possibly  occurred  in  active  subjects  as  manifested  by 
increased  interstitial  extra-axonal  water  and  delayed  clear¬ 
ance  of  Gd-DTPA.  These  results  would  be  consistent  with 
a  decrease  in  diffusivity  with  histopathological  evidence  of 
axonal  injury  in  rats  following  traumatic  brain  injury.37,38 


Average  1st  and  2nd  Imaging 


Time(s) 

— Avg  1st  Imaging  Observed  #  Avg  1st  Imaging  Function  — +—  Avg  2nd  Imaging  Observed  — Avg  2nd  Imaging  Function 

FIGURE  2.  Comparison  of  actual  observed  group  data  to  exponential  function  for  both  the  first  and  second  MRIs.  On  both  occasions,  the  actual  data 
tightly  follow  the  derived  exponential  function  (%2,  p  =  0.2395;  Spearman’s  correlation,/?  <  e-16),  permitting  cross  group  comparison  of  the  regression  equa¬ 
tion  parameters. 
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TABLE  III.  Nonlinear  Regression  Equation  Parameters  and  Washout  Rate  and  Ratio 


Group 

MRI  1  (Pre-exposure) 

MRI  2  (Postexposure) 

t  Value//?  Value 

A 

B 

A 

B 

A 

B 

Sham 

0.087  ±  0.002 

0.439  ±  0.030 

0.085  ±  0.016 

0.435  ±  0.078 

0.217/0.86 

0.095/0.91 

Study 

0.100  ±0.014 

0.342  ±0.113 

0.098  ±  0.008 

0.336  ±  0.056 

0.269/0.80 

0.107/0.93 

Washout  Rate 

Intensity  Ratio 

Washout  Rate 

Intensity  Ratio 

Washout  Rate 

Intensity  Ratio 

Sham 

0.262  ±  0.048 

0.995  ±  0.041 

0.265  ±  0.090 

0.975  ±  0.032 

-0.047/0.97 

0.679/0.68 

Study 

0.262  ±  0.048 

0.262  ±  0.048 

0.262  ±  0.048 

0.262  ±  0.048 

0.180/0.79 

0.895/0.34 

A  and  B  are  fit  coefficients  from  Intensity  =  A  x  ln(Y)  +  B  where  t  =  time.  Intensity  ratio  =  final/base  intensity. 


Repeated  anesthesia  and  MRI  in  the  sham  animals  did  not 
produce  any  significant  changes  on  these  indices  sensitive 
to  WM  integrity.  This  suggests  the  observed  MRI  changes 
were  possibly  related  to  the  hypobaric  exposure,  although 
other  contributing  factors  cannot  be  excluded. 

The  lack  of  new  WMH  lesions  in  this  model  is  consis¬ 
tent  with  similar  findings  of  an  in-progress  single  exposure 
human  study.  We  posit  that  the  initial  injury  is  a  diffuse 
process  and  that  subsequent  FLAIR  change  occurs  in  those 
instances  where  repair  has  been  inadequate.  With  this 
hypothesis,  initial  injury  might  be  reflected  by  BBB  integ¬ 
rity  and  interstitial  fluid  content  changes,  although  FLAIR 
change  would  not  be  expected.  Although  the  Mu  trend  data 
and  its  relationship  to  Gd-DTPA  washout  rate  are  encourag¬ 
ing,  the  lack  of  significant  difference  in  other  parameters  sug¬ 
gests  this  pig  model  of  cerebral  injury  from  hypobaria  needs 
further  development. 

One  substantial  difference  between  this  model  and  occu¬ 
pational  hypobaric  exposure  in  humans  is  that  the  animals 
were  exposed  while  under  anesthesia.  Humans  are  awake 
and  active  during  occupational  exposure,  and  physical 
activity  at  altitude  is  known  to  potentiate  brain  injury.39  In 
addition,  anesthesia  may  act  as  a  neuroprotective  agent, 
thus  making  the  animal  experiments  less  realistic.  Necropsy 
suggested  the  high  mortality  rate  in  this  study  was  second¬ 
ary  to  the  combination  of  decompression  sickness  and/or 
ventilator-associated  barotrauma  with  severe  pulmonary 
injury  and  sepsis.  Although  neuropathological  abnormali¬ 
ties  were  observed,  the  relationship  to  possible  hypobaric- 
induced  injury  is  unclear.  These  issues  are  being  addressed 
by  utilizing  a  nonanesthetized,  physically  active  pig  in  a 
phase  II  trial.  Additionally,  this  phase  II  model  utilizes 
an  expanded  set  of  quantitative  structural  (T1  and  T2 
weighted)  and  physiological  (micro structural  properties  of 
molecular  diffusion,  DSC,  and  concentrations  of  neuro¬ 
chemicals)  parameters  for  a  more  comprehensive  assess¬ 
ment  of  the  pathophysiological  process. 

The  consistency  of  measurements  in  this  model  suggests 
that  MRI  can  be  used  in  a  miniature  swine  model  to  study 
brain  trauma  induced  not  only  by  hypobaric  exposure  but 
resulting  from  other  factors.  Susceptibility- weighted  perfu¬ 
sion  MRI  may  be  difficult  to  obtain  in  the  acute  TBI  setting 
secondary  to  patient  contraindications  and  is  limited  by 


signal  loss  in  the  presence  of  blood  products,  but  may  be 
valuable  in  subacute  or  chronic  TBI  particularly  when  com¬ 
bined  with  other  conventional  and  advanced  MRI  techniques 
including  DTI.40-42  Patients  with  mild  TBI  and  normal  con¬ 
ventional  imaging  have  demonstrated  perfusion  deficits  that 
correlate  with  neurocognitive  impairment.43,44  TBI  causes 
axonal  and  capillary  disruption45  that  can  be  observed  as 
WMH  on  FLAIR  MRI.44,46-48  DTI  can  detect  microscopic 
tissue  damage  in  TBI  and  other  central  nervous  system  dis¬ 
orders,  particularly  those  that  affect  WM,  and  provide  markers 
that  may  be  capable  of  predicting  clinical  outcome.40,41,43,44 
Thus,  although  the  focus  of  this  study  has  been  on  brain  injury 
following  hypobaric  exposure,  the  findings  from  this  model 
are  likely  to  be  relevant  for  mild-to-moderate  TBI,  thus 
avoiding  the  focal  cerebral  injury  associated  with  other  ani¬ 
mal  models.  This  model  suggests  a  general  disruption  of 
axons  has  occurred,  possibly  similar  to  the  shear  injury  in 
TBI.  Thus,  potentially  this  model  might  provide  insight 
into  the  injury  biomechanics  and  axonal  pathophysiology 
of  TBI,  hence  addressing  one  of  the  major  contributors  to 
subsequent  disability. 
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